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Abstract 

The transonic performance of a dual-throat, single- 
expansion-ramp nozzle (SERN) was investigated with 
the PARC computational fluid dynamics (CFD) code, 
an external flow Navier-Stokes solver. The nozzle 
configuration was from a conceptual Mach 5 cruise 
aircraft powered by four air-breathing turboramjets. 
Initial test cases used the two-dimensional version of 
PARC in Euler mode to investigate the effect of 
geometric variation on transonic performance. 
Additional cases used the two-dimensional version in 
viscous mode and the three-dimensional version in both 
Euler and viscous modes. Results of the analysis 
indicate low nozzle performance and a highly three- 
dimensional nozzle flow at transonic conditions. In 
another comparative study using the PARC code, a 
single-throat SERN configuration for which 
experimental data were available at transonic 
conditions was used to validate the results of the 
over/under turboramjet nozzle. 

Symbols 

Cf force coefficient normalized by the ideal 
force 

P 0 ambient static pressure 
P s ramp static pressure 
X axial distance on ramp surface 
Y normal distance on ramp surface 

Subscripts: 

g gross 
x axial 
y normal 
int internal 
ramp nozzle ramp surface 
cowl cowl surface 

* 

Nozzle Technology Branch. 

This paper is declared a work of the U.S. Government 
and is not subject to copyright protection in the United 
States. 


Introduction 

Recent advances in technologies required for 
hypersonic vehicles indicate that Mach 4 to 6 vehicles 
are now feasible. These advances have renewed interest 
in the development of hypersonic vehicles for both 
military and commercial applications. A new program 
was developed recently from the Inlets and Nozzle 
Concepts for Advanced Air-breathing Propulsion 
Systems study sponsored by the Air Force. This new 
joint program was conducted by the NASA Langley 
and NASA Lewis Research Centers to evaluate a 
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Mach 5 waverider cruise vehicle (shown in Figure 1). 
Langley conducted the vehicle design and mission 
analysis activities, whereas Lewis estimated the 
installed propulsion system performance. One of the 
objectives of this program was to develop performance 
prediction methods for single- and dual-flow exhaust 
nozzle systems. 

The propulsion system for the Mach 5 waverider 
used four of the over/under airbreathing turboramjet 
engines shown in Figure 2. For these types of 
propulsion systems, the exhaust nozzles must operate 
over nozzle pressure ratios (NPR, internal total 
pressure divided by ambient static pressure) varying 
from about 2 to 300. A complexity arising from this 
propulsion concept is that both the turbojet and ramjet 
operate together, producing dissimilar nozzle flows, 
over a wide range of Mach numbers. 

The wide operating range of this system requires a 
variable throat in each nozzle. Flaps or other 
mechanical devices can be used to vary the nozzle area, 
but they add weight to the exhaust system. One of the 
attractive advantages of using a single-expansion-ramp 
nozzle (SERN) instead of a conventional, two- 
dimensional, convergent-divergent nozzle is that the 
underside aft region of the aircraft can be used as the 
expansion surface for the flow. This helps to reduce the 
overall weight of the vehicle. However, previous 
studies have shown that hypersonic SERN nozzles of 
this type generally have low thrust performance at 
transonic and low supersonic Mach numbers because 
the nozzles are overexpanded. 
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The objective of the current work was to obtain 
initial nozzle performance data for the Mach 5 
waverider and to optimize the nozzle at transonic 
conditions. The following sections present the 
procedures used and results obtained from this analysis. 

Nozzle Geometry and Flow Conditions 

Figure 3 shows the dual-throat turboramjet nozzle 
geometry used for this study, with its assumed variable 
geometry features. The RAO code (method of 
characteristics) was used to design and optimize the 
isentropic ramp contour at a design point of Mach 5 
with the ramjet only. On the basis of this design point, 
the nozzle had a ramp length of 465 in. and an overall 
area ratio of 15. The flap position separating the two 
internal flows at Mach 1.2 was based on the ramp 
contour at a design point of Mach 5 and the required 
throat areas at Mach 1.2 from the engine cycle deck. 
The throat heights for the turbojet and ramjet were 
determined to be 10.3 and 19.4 in., respectively, and 
the analysis assumed a rectangular cross-sectional 
throat area with a total width of 180 in. (45-in. wide for 
each of the four engines). The turbojet had an internal 
area ratio from the throat to the end of the splitter of 
1.92, whereas the ramjet had a ratio of 1.75. 


Table 2 Coordinates for bottom 
surface of turbojet 
[All dimensions given in feet.] 


X 

Y 

3.250 

3.042 

3.432 

3.031 

3.613 

3.020 

3.794 

3.010 

3.975 

3.000 

4.156 

2.989 

4.337 

2.979 

4.518 

2.969 

4.698 

2.958 


X 

Y 

4.879 

2.948 

5.058 

2.937 

5.237 

2.927 

5.417 

2.917 

5.517 

2.902 

5.617 

2.887 

5.717 

2.872 

5.817 

2.857 

5.917 

2.842 


Table 3 Coordinates for upper surface of ramjet 
[All dimensions given in feet] 


X 

Y 

-0.333 

1.617 

.055 

1.757 

.438 

1.891 

.812 

2.010 

1.170 

2.118 

1.511 

2.212 

1.831 

2.316 

2.127 

2.401 

2.400 

2.466 

2.648 

2.506 


X 

Y 

2.872 

2.541 

3.072 

2.570 

3.250 

2.596 

3.432 

2.622 

3.613 

2.647 

3.794 

2.671 

3.975 

2.691 

4.156 

2.710 

4.337 

2.727 

4.518 

2.741 


X 

Y 

4.698 

2.754 

4.879 

2.767 

5.058 

2.778 

5.237 

2.791 

5.417 

2.802 

5.517 

2.809 

5.617 

2.815 

5.717 

2.821 

5.817 

2.827 

5.917 

2.833 


Tables 1 to 3 give the overall coordinates for the 
nozzle geometry. The inflow total and freestream 
conditions were based on engine data obtained from the 
cycle deck program. For the turbojet, the total pressure 

and temperature were 6778 lb/ft 2 and 3897 °R, 
respectively. For the ramjet, the total pressure and 
temperature were 1334 lb/ft 2 and 3835 °R, 
respectively. The freestream static pressure and 
temperature were 628 lb/ft" and 412 °R, respectively. 
The turbojet operated at an NPR of 10.8, whereas the 
ramjet operated at an NPR of 2. 1. 


Table 1 Coordinates for upper surface of turbojet 
and expansion ramp 
[All dimensions given in feet] 


X 

Y 

26.946 

10.136 

27.929 

10.279 

28.936 

10.418 

29.964 

10.553 

31.012 

10.683 

33.158 

10.929 

35.356 

11.152 

36.469 

11.255 

37.588 

11.352 

38.709 

11.442 


X 

Y 

13.572 

7.315 

14.505 

7.580 

15.540 

7.858 

16.683 

8.149 

17.940 

8.449 

19.313 

8.757 

20.806 

9.068 

22.418 

9.380 

23.268 

9.535 

25.053 

9.840 


X 

Y 

3.250 

3.900 

4.337 

4.052 

5.417 

4.250 

6.546 

4.779 

7.536 

5.213 

8.637 

5.662 

9.736 

6.062 

10.458 

6.315 

11.320 

6.615 

12.349 

6.951 


Computational Methods 
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For the present study, the PARC code was used on 
the CRAY-YMP to determine the nozzle performance 
at Mach 1.2. PARC, an external flow Navier-Stokes 
solver based on the ARC code, can analyze a variety 
of propulsion flows. It can be operated in either Euler 
mode or full Navier-Stokes mode. PARC also allows 
for a variety of boundary conditions, making it a very 
versatile flow solver. For this study, the Thomas 
algebraic model was used as the turbulence model. 

Figure 4 shows a typical two-dimensional 
computational grid generated by 3D INGRID and 
used for this study, and Figure 5 shows a typical 
three-dimensional grid. The two-dimensional grids 
had about 18 000 grid points, whereas the three- 
dimensional grids had about 500 000 grid points. Only 
the grid points of the supersonic portion of the nozzle 
from the throat to exit were used in the computation. 
The analysis neglected the inner sidewalls separating 
the four engines for the three-dimensional cases. 
Upstream inflow conditions began at the throats of 
each nozzle and were fixed at their appropriate 
pressures and temperatures, whereas the downstream 
outflow condition was extrapolated. All the surfaces 
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were considered to be slip walls for the Euler cases and 
to be no-slip, adiabatic walls for the viscous cases. The 
free stream boundary condition was used for the 
external flow of Mach 1.2. 

A post-processing program calculated the 
parameters used to quantify the nozzle 
performance — the axial and normal force coefficients, 

and c f.r Figure 6 shows the control volume used 
for the analysis. For each directional force, the program 
calculated the internal stream thrust by integrating the 
momentum and pressures at the common exit plane. It 
calculated the ramp pressure force in a similar way by 
integrating the pressures on the ramp surface. The cowl 
boat-tail force was estimated assuming a Prandlt-Meyer 
expansion. These three individual force components 
(internal stream thrust, ramp force, and cowl force) all 
contributed to the overall gross force. The ideal force 
for each nozzle flow is a function of the throat area, 
NPR, and specific heat ratio. For the specific heat ratio, 
the average specific heat ratios of the turbojet and the 
ramjet were used. The overall ideal force was the sum 
of the individual ideal forces for the two nozzles. When 
the ideal force calculations changed less than 1 percent, 
the solution was considered to be converged. On the 
average, it took about 20 000 iterations for each 
solution to converge. 

Results 

The first configuration analyzed with the PARC2D 
code in Euler mode was the baseline case, where the 
cowl angle was set at 0°. Figure 7 shows a pressure 
mismatch between the exit planes of the turbojet and 
the ramjet. The turbojet flow is the dominant flow in 
this dual mode operation during transonic flight. As 
mentioned before, the turbojet is operating at an NPR 
of around 1 1 , compared with an NPR of 2 for the 
ramjet Because of the low ramjet exit pressure, the 
flow from the turbojet tends to expand locally into the 
ramjet flow. On the other hand, the ramjet flow is 
overexpanded, and an oblique shock appears at the 
trailing edge of the flap, splitting the two flows. Also, 
there seems to be a small local normal shock causing a 
subsonic region in the flow near the cowl exit of the 
ramjet 

In the axial direction (at a 0° cowl angle), the 
gross force coefficient calculated was about 
0.722, indicating low nozzle performance. However, 
the internal stream thrust coefficient C^ mt was quite 
high at 0.934. The penalty came from the external 


expansion ramp, which produced a negative ramp force 
coefficient C 'f tXjT3mp of -0.212 because of the 
overexpansion of the nozzle at Mach 1.2. In the normal 
direction, the gross force coefficient Cf y g was about 
-0.847. Similar effects were seen in this component as 
well where the expansion ramp produced a large 
negative force: Cj y >ramp was about -0.829, whereas 
Cfy/mt was on ty about -0.018. The cowl boat-tail force 
was assumed to be zero in both directions since the 
cowl had 0° deflection. 

A first approach to improve the nozzle 
performance was to reduce the pressure difference at 
the exits of the two flows by increasing the local 
pressure that the turbojet was expanding into. The cowl 
of the ramjet was deflected upwards 5° to increase the 
ramjet exit pressure by decreasing the exit area. A 
normal shock appeared inside the ramjet which helped 
mi nimi ze the pressure difference between the two 
flows due to the increase in ramjet exit pressure 
(Fig. 8). Consequently, C fx.g improved slightly from 
the baseline configuration but was still low at 0.755, 
and Cf ^ mi increased from 0.934 to 0.953. The cowl 
deflection slowed down the expansion of the turbojet 
flow on the ramp, thus lowering C to about 
-0.188 from -0.212. However, the cowl deflection also 
caused a negative boat-tail force with a coefficient of 
-0.010. As a result, the net effect in the axial force was 
a small increase over the baseline configuration. The 5° 
cowl deflection had a greater effect in the individual 
forces in the normal direction. The normal internal 
stream thrust coefficient, Cf y inV increased from -0.018 
to -0.121, and Cf y T3mp decreased from -0.829 to -0.771. 
Also C f y cowl increased from 0.000 to -0.1 12. Overall, 
the gross normal force Cj y g increased from -0.847 to 
-1.004. 

Another configuration was developed by 
deflecting the cowl at an angle that would equalize the 
pressures exiting from both nozzles. A one-dimensional 
analysis was used to determine the exit area from the 
exit static pressure (which equals the exit static 
pressure of the turbojet), total pressure, and throat area 
of the ramjet The cowl angle associated with this new 
ramjet exit area was determined to be deflected 
upwards about 8.4°. The oblique shock at the end of 
the flap splitting the two flows was eliminated, and the 
turbojet flow no longer expanded into the ramjet flow 
as shown in Figure 9. The flow in the ramjet, however, 
stayed subsonic for the most part and did not become 
supersonic until it interacted with the freestream flow 
at the cowl exit. 
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Once again, no significant changes were seen in 
the axial direction. The internal stream thrust, ramp 
pressure force, and boat-tail force coefficients were 
0.961, -0.177, and -0.026, respectively. These three 
coefficients contributed to a C fx.g of 0.758, which was 
a slight increase from the baseline configuration at 
0.722 but no dramatic change from the 5° deflection at 
0.755. The changes were more apparent for the normal 
forces: Cf yin{ increased to -0.153, Cf y T2mp decreased 
to -0.736, and C f y ;cowI increased to A). 174. However, 
these changes were not as significant when compared 
with the 5° cowl deflection. Overall, the gross normal 
force coefficient C f yg was about -1.063. 

Figure 10 compares the pressure distributions 
P s /P 0 (ramp static pressure divided by ambient static 
pressure) along the expansion ramp at the cowl angle 
settings given in Figures 7 to 9 and Tables 4 and 5. The 
static pressures along the ramp did increase slightly as 
the cowl angle increased. This increase in pressure was 
reflected in the increase in the force coefficients. The 
effect was more significant in the normal direction than 
in the axial direction. However, the static pressures 
along the expansion ramp still stayed below the 
ambient static pressure. Tables 4 and 5 summarize the 
two-dimensional Euler performance calculations at the 
various cowl angles. 

Because the ramp pressures always remained 
below the ambient pressure in the two-dimensional 
analysis, the next step was to investigate the 
three-dimensional relieving effect for the baseline 
configuration with PARC3D in Euler mode. The 
modeling of the freestream flow came into the 
calculation just at the cowl exit and downstream of 
that; this minimized the computational time. Also, only 
half the nozzle was modeled because the nozzle was 
symmetric along the expansion ramp. PARC3D 
predicted a significant pressure rise along the centerline 
of the expansion ramp, as shown in Figure 11, whereas 
the two-dimensional case always stayed below the 
ambient static pressure. Figures 12 and 13 show the 
Mach number contours along the centerplane and on 
the surface of the expansion ramp. The pressure rise 
tends to decrease away from the centerline. At 
transonic speed, the freestream helped pressurize the 
flowfield because the nozzle was highly overexpanded. 
The pressurization is reflected in the ramp force 
calculations. 

The PARC3D configuration in Euler mode had a 
Cfx,g of 0.842 — significantly higher than the 
two-dimensional prediction of 0.722 — and had 


Table 4 Axial force coefficients for 
two-dimensional Euler analyses 


Cowl 

angle, 

deg 

Axial force coefficients 

CfixjM 


Cfo oowl 

Cfx,g 

0.0 

0.934 

-0.212 

0.000 

0.722 

5.0 

.953 

-.188 

-.010 

.755 

8.4 

.961 

-.177 

-.026 

.758 


Table 5 Normal force coefficients for 
two-dimensional Euler analyses 


Cowl 

angle, 

deg 

Normal force coefficients 

Cfy.'m 

Cfy.nmp 

f'f.y.mwl 

Cf.y.g 

0.0 

-0.018 

-0.829 

0.000 

-0.847 

5.0 

-.121 

-.771 

-.112 

-1.004 

8.4 

-.153 

-.736 

-.174 

-1.063 


individual axial components, and of 

0.951 and -0.109, respectively. The performance in the 
normal direction was greatly improved: Cj y mi was 
calculated to be -0.027, and Cf y TZJttp was -6.241, with 
an overall c /» g of -0.268. Once again, the cowl boat- 
tail force was assumed to be zero in both directions for 
the 0° cowl deflection. The axial and normal internal 
stream thrust coefficients, Cy^ mt and C f y inv varied 
slightly between the two-and three-dimensional cases 
because different grid densities were used to minimize 
computation time. Ideally, the values would be the 
same for both cases. 

The next approach in the study was to look at the 
viscous effects for the baseline configuration with 
PARC in both two and three dimensions. The grids 
were modified for the viscous configurations with 
tighter packing in the shear layer and near the wall 
surfaces. Figure 14 shows the Mach number contours 
for the two-dimensional viscous case. A boundary layer 
can be seen along the external expansion ramp. The 
flow seems to be separating near the exit of the ramjet 
both at the top and bottom surfaces. Also, the flow 
from the turbojet does not expand into the ramjet as 
much as in the two-dimensional Euler case at a cowl 
angle of 0°. For the force calculations in the axial 
direction, Cf ^ mi was about 0.947 and was 

about -0.205, to give a C^ g of 0.742. In the normal 
direction, Cf yiDt and C fy T3mp were determined to be 
-0.062 and -6.825, respectively. The net force 
coefficient Cj y g was -0.887. The force coefficients 
were, however, slightly higher than for the 
two-dimensional Euler case. One reason might be that 
the grid was not packed closely enough near the wall 
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surfaces to accurately capture the boundary layer. 
Overall, there were no real significant performance 
differences in the forces in comparison to the 
two-dimensional Euler case. 

Figure 15 shows the centerplane Mach number 
contours for the three-dimensional viscous case. This 
figure is very similar to the two-dimensional viscous 
case with the exception that the flow separates along 
the expansion ramp. Figure 16 shows the pressure 
distribution comparison for the baseline configuration. 
The pressure rise for the three-dimensional viscous 
case was significantly higher than for the 
two-dimensional cases but slightly lower than for the 
three-dimensional Euler case. For the performance 
calculations in the axial component, C u int was °- 906 
and Cf ^ ramp was -0.086, to give a Cf Ag of 0.820. In the 
normal component, Cy y in( was -0.078 and C^ y T3mp was 
-0.187, to give a g of -0.265. In both components, 
C/cowi was assumed to be zero. The three-dimensional 
calculations were much more accurate than the two- 
dimensional calculations, but there were no major 
differences between the Euler and viscous calculations. 
Tables 6 and 7 summarize the comparison between the 
Euler and viscous calculations for the baseline 
configuration. 

Evaluation of Methodology 

There was an uncertainty in the significant 
pressure rise for the over/under nozzle concept 
predicted by the three-dimensional PARC analysis. 
Another comparative study was done with a similar 
type of hypersonic nozzle for which experimental data 


Table 6 Axial force coefficients for baseline configuration 

[C w = 0] 



Mode 

Axial force coefficients 

Cf.xjni 

^Uramp 

C Ug 

PARC2D 

Euler 

0.934 

-0.212 

0.722 

PARC2D 

Viscous 

.947 

-.205 

.742 

PARC3D 

Euler 

.951 

-.109 

.842 

PARC3D 

Viscous 

.906 

-.086 

.820 


Table 7 Normal force coefficients for baseline configuration 
[C/y.cowl ~ Qj 



Mode 

Normal force coefficients 

C fyM 

^fy.nm p 

Cf.y.g 

PARC2D 

Euler 

-0.018 

-0.829 

-0.847 

PARC2D 

Viscous 

-.062 

-.825 

-.887 

PARC3D 

Euler 

-mi 

-.241 

-.268 

PARC3D 

Viscous 

-.078 

-.187 

-.265 


were available at transonic conditions to validate the 
PARC results. The turbulence model used for this study 
was, once again, the Thomas algebraic model. The 
experimental work in reference 8 was conducted in the 
8- by 6-Foot Supersonic Wind Tunnel at NASA Lewis 
on a National Aerospace Plane (NASP) type SERN 
with a single flow path (Figure 17). The nozzle, which 
was enclosed by sidewalls, was 8-in. wide with a 1-in.- 
high throat that expanded to an internal area ratio of 
1.29. The external expansion ramp was 8-in. wide and 
had an initial angle of 17° and a trailing edge angle of 
8°. The experimental conditions included various 
tunnel Mach numbers and NPR's, cold and hot flows, 
with and without external burning, and different nozzle 
configurations. The data point chosen was at a tunnel 
Mach number of 1.2 and an NPR of 10 for cold flow 
without external burning. 

Experimental data showed a significant 
three-dimensional effect along the centerline of the 
expansion ramp at Mach 1.2 where the freestream 
interacted with the primary flow (Figure 18). The 
pressure rise decreased away from the centerline 
similar to the over/under turboramjet configuration. 
Once again, both the two- and three-dimensional 
versions of PARC were used in the computation. The 
PARC2D results underpredicted the static pressures 
along the expansion ramp. The three-dimensional 
results matched fairly well with the experimental data 
for pressures along the centerline. However, away from 
the centerline, PARC3D's predictions did not agree 
well with the experimental data. Also the overall 
pressure rise in the flowfield was not as high as the 
waverider’s results, but nevertheless a 
three-dimensional effect was visible. 

Summary and Conclusions 

Several configurations of an over/under 
turboramjet nozzle concept were investigated with the 
PARC code in Euler mode and viscous mode with the 
Thomas model. Computational fluid dynamics (CFD) 
results show low-thrust performance at transonic flight 
conditions and a highly three-dimensional flowfield. 

The three-dimensional results agree qualitatively with 
the experimental data for the hypersonic nozzle at 
Mach 1.2. PARC accurately predicted the pressure rise 
along the centerline of the nozzle. The cowl deflection 
had a minimal benefit in the axial direction but had a 
larger negative effect in the normal direction. There 
were no significant differences in the nozzle 
performance calculations of the Euler and the viscous 
cases. Future CFD studies on these types of nozzles 
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should utilize a three-dimensional analysis to 
accurately capture the three-dimensional effects at 
transonic conditions. An Euler analysis is sufficient for 
an initial parametric study. 
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Figure 1 . — Three views of the Mach 5 waverider cruise vehicle. 2 (All dimensions given in inches unless 
marked otherwise.) 
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Figure 2. — Over/under turboramjet propulsion system. 
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Figure 3. — Dual-throat turboramjet nozzle geometry. 
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Figure 4. — Typical two-dimensional computational grid used for this study. 
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Figure 5. — Typical three-dimensional computational grid used for this study. 
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Figure 6. — Control volume used for force calculations. 
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Figure 7. — Mach number contours for two-dimensional case in Euler mode at 0° cowl angle. 



Figure 8. — Mach number contours for two-dimensional case in Euler mode at 5° cowl angle. 
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Figure 9. — Mach number contours for two-dimensional case in Euler mode at 8.4° cowl angle. 
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Figure 1 0. — Comparison of ramp pressure distributions for two-dimensional 
case in Euler mode at various cowl angles. 



Figure 1 1 . — Comparison of ramp pressure distributions for two- and three- 
dimensional cases in Euler mode. 
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Figure 12. — Mach number contours along centerplane for three-dimensional case in Euler mode. 
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Figure 13. — Mach number contours along expansion for three-dimensional case in Euler mode. 
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Figure 14. — Mach number contours for two-dimensional case in viscous mode. 


Mach number 
contour 


3.2 


0.0 




Figure 15. — Mach number contours along centerplane for three-dimensional case in viscous mode. 
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Figure 16. — Comparison of ramp pressure distributions for baseline configuration. 
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Figure 17. — Single-flow Single-expansion-ramp nozzle (SERN) configuration 
in the 8-Foot by 6-Foot Supersonic Wind Tunnel. 
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Figure 1 8. — Comparison of ramp pressure distributions for National Aerospace 
Plane (NASP) configuration at tunnel Mach number of 1 .2 and NPR of 1 0. 
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